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Abstract. 

We present the discovery of nine quasars behind the Large Magellanic Cloud, with emission redshifts ranging from 0.07 to 2.0. 
Six of them were identified as part of the systematic variability-based search for QSOs in the objects from the OGLE-II database. 
Combination of variability-based selection of candidates with the candidates' colours appears to be a powerful technique for 
identifying quasars, potentially reaching ~50% efficiency. We report an apparent correlation between variability magnitude 
and variability timescale, which - if confirmed - could put even more constraints on QSO candidate selection. The remaining 
three quasars were identified via followup spectroscopy of optical counterparts to X-ray sources found serendipitously by the 
Chandra X-ray Observatory satellite. Even though the locations of the candidates were quite uniformly distributed over the 
LMC bar, the confirmed QSOs all appear near the bar's outskirts. 
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1. Introduction 

Recent advances in observational techniques and the availabil- 
ity of large variability monitoring databases have opened up a 
possibility for searches for quasars behind dense stellar fields. 
Until recently, those areas were avoided in QSO surveys, pri- 
marily because of candidate confusion problems. At the same 
time, such quasars are very intere sting, e.g. a s background 
sources for absorption studies (e.g. lProchaska et alJl2002l). as 
fixed reference points for dynamical studies ("e.g. lAneuita et al] 
2000), or even as means to determine very acc urate distances 
to the foreground objects (foraine & Bondl2004h . The launch of 
the Chandra X-ray Observatory - with its superb pointing ac- 
curacy and spatial resolution - enabled X-ray based searches, 
while microlensing campaigns (such as MACHO and OGLE) 
provided the possibility for variability-based identifications. 

The fields of particular interest are nearby galaxies, such as 
the Magellanic Clouds. Application of new techniques resulted 
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in the number of known quasars in the general direction of the 
Clouds more than trebling since 2002 JPobrzvcki etafll2002l 
2003 JS lGeha et all2003l) . and this difference is more dramatic 
in the direction of the inner parts of the Clouds, where almost 
all presently known QSOs were discovered in recent years. 

Despite the remarkable success shown by these methods, 
there remained an area which was still mostly unexplored: the 
bar of the Large Magellanic Cloud, where only four quasars 
were known. At the same time, this was the region where find- 
ing QSOs would be most interesting for studying the dynam- 
ics of the LMC. The expected proper motion at the distance to 
the LMC is of the o rder of only few miliarcseconds per year 
jAnguita et alJ200ol and references therein). To study bulk mo- 
tions it is thus critical to have a sizable number of fixed refer- 
ence points placed in "strategic" positions, preferably behind 
the bar. 

There were two main reasons for the scarcity of known 
QSOs in this region. First, the X-ray-based searches - which 
are responsible for three out of four quasars behind the LMC 
bar - are limited by the relatively small number of Chandra 
pointings towards the LMC; on top of that the known fields 
have been only partially explored. Second, the only variability- 
based sample that has so far been followed up spectroscopi- 
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Table 1. New quasars. 



J2000.0 coordinates 


Zera 


V [mag] 


Sel. 


ID" 


J0505 14.37-700050.9 


1.83 


20.08 


Var. 


L98 


J05 1244.64-693857.2 


0.31 


19.51 


Var. 


L81 


J051433.58-693245.7 


1.37 


18.87 


Var. 


L84 


J052141. 75-703029.2 


1.17 


19.31 


Var. 


L154 


J052627.69-693424.7 


2.09 


19.81 


Var. 


L24 


J053400.06-703 127.9 


1.14 


19.39 


Var. 


LI 


J050550.35-675017.5 


0.07 


19.00 


X-ray 




J050900.76-672114.1 


1.63 


20.58 


X-ray 




J050940.60-685320.6 


0.67 


19.79 


X-ray 





" Candidate identification in lEveJ2002l 



cally is the MACHO sample Jdeha et alJl2003l) . This sample - 
although it is the biggest contributor to the number of known 
quasars in the vicinity of the LMC - has not been followed up 
completely, and the completed part is somewhat biased towards 
the outer parts of the LMC. 

In this paper we present the results of the search for QSOs 
behind the Large Magellanic Cloud. The variability-based 
identi fication of c andidates in the OGLE-II dUdalski et alJ 
2000; IZebrun et al]l200ll) data, which cover 4.6 deg 2 co ntain- 
ing the bar of the LMC, was performed by Ever (2002), and 
this list was the main part of our candidate pool. We also fol- 
lowed up on several can didates from the X-ray selected sample 
JPobrzvcki et alJ2002h . 

This work is a part of an ongoing project ( Dob rzvcki et alJ 
20021 l2003albl) . The variability-based technique in the vicinity 
of the Small Magellanic Cloud turned out to be very efficient, 
nearing 40%, although the sample is rather small (five quasars 
among seventeen candidates). Due to several effects - some re- 
lated to the intrinsic differences between the Clouds and some 
related to the nature of the OGLE data - the list of candidates 
in the direction of the LMC was expected to be more contam- 
inated with non-quasars. See Sections 2 and 3 in lEverl J2002I) 
for a discussion of those effects. 



2. New quasars 

2.1. Observations and identifications 

The spectroscopic observations of the candidates were per- 
formed on several nights using two instruments. On 2002 
January 22-23, 2002 November 07 and 17 we used the 
Magellan Baade 6.5-meter telescope, the LDSS-2 imaging 
spectrograph with 2048x2048 SITe#l CCD, 1.03 arcsec slit 
and 300 1/mm MedBlue grism, yielding a nominal resolution of 
13.3 A. On 2004 December 13-17 we used the ESO 3.6-meter 
telescope, the EFOSC2 imaging spectrograph and Loral/Lesser 
2048x2048 CCD#40 detector, 1 arcsec slit and 300 1/mm grism 
#06, yielding a nominal resolution of 12.9 A. 

During all nights, we observed Ham uv et alJ i 19921) spec- 
trophotometric standards and took exposures of He-Ne-Ar 
lamps for wavelength calibrations. All spectra were reduced in 
the standard way using IRAF. 



The followup observations revealed nine previously un- 
known QSOs. Their rest frame spectra are shown on Figure^ 
and a summary of their properties is shown in Tabled 

2.2. Variability selected quasars 

The set of candidates was constructed f rom the objects in 
the OGLE-II database of variable objects (Zeb ruri et alJl200ll) 
using several criteria. In addition to colour-magnitude and 
colour-colour criteria, which removed many variable stars with 
known characteristics (e.g. Cepheids), a more elaborate crite- 
rion, aimed at selecting objects which (like QSOs) show larger 
variability on longer time scales, was applied. This was done by 
calculating variogram, a function showing time lag versus the 
spread of all magnitude differences than the lag. The variogram 
is steeper for objects which behave in the desired way, and an 
empirical cutoff o n the v ariogram slope for long timescales was 
applied. See lEverlll2002l) for details. 

In total, the list of variability-selected QSO ca ndida t es con - 
tained 118 objects. One of them, object L92 in lEveJ j2002t) . 
was identified as a z e m — -61 QSO in the X-ray-based search 
of iDobrzvcki et alJ J2002I) . Of remaining objects, six were 
als o identified as Q SO candidates in the MACHO database 
bv lGehaetall(l2003l) . Two of those objects, LI 14 and L155, 
were identified as quasars at, respectively, z em = 1.81 and 
Zem = 0.90; the other four objects turned out to be stars. 

During our observing runs, we completed observations of 
108 objects out of the remaining 111 candidates, identifying 
six of them as QSOs. OGLE light curves of those six quasars 
are shown in Figure|2 We do not see qualitative differences be- 
tween quasars and stars in the appearances of their light curves, 
but see Section l3~3*l below for discussion of possible quantita- 
tive difference between them in how variogram slope relates to 
variability scale. 

We consider the sample to be fully observed, since the re- 
maining three objects are very unlikely to be quasars, given 
their distinctly non-QSO colour properties (see Section l3~TV 

Combined with three quasars identified previously, we 
found the original variability-selected sample to include nine 
identified QSOs, i.e. a ~8% success rate. This is lower that 
the rate seen in a simil ar search in the dire ction of the Small 
Magellanic Cloud dDobrzvcki et al.ll2003al) but, as mentioned 
earlier, not unexpected. The variability-selected quasars have 
emission redshifts ranging from 0.31 to 2.09, and their discov- 
ery spectra (Fig.^ seem to show only typical quasar features. 

2.3. X-ray selected quasars 

During the runs we took sev eral spectra of the X-ra y selected 
candidates from the list from Dobrzvc ki et al.N2002l) . We refer 
the reader to this paper for detailed description of the sample; 
we only mention here that this sample is in no way homoge- 
neous, therefore we will not derive any quantitative properties 
regarding the efficiency of the method or completeness of the 
resulting quasar harvest. 

In total, we observed 30 objects. Of those, three turned out 
to be quasars, at emission redshifts of 0.07, 0.67 and 1.63. We 
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Fig. 1. Rest-frame discovery spectra of quasars presented in this paper. First two rows show variability-selected QSOs, third row 
shows X-ray-selected objects. 
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Fig. 2. OGLE (Zeb run etall200ll) li ght curves for six variability-selected quasars (top six in Figure^. MJD 50,400 corresponds 
to UT 1996 November 13. 
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note here that the first of them, QSO J050550.35-675017.5, is 
a relatively strong and highly variable X-ray source (observed 
Chandra ACIS-7 count rate of 0.12-0.38 counts/sec), and per- 
haps could be used for geometrical LMC distance determina- 
tion using the method described in lDraine & Bond J2004I) . al- 
though it is likely that a next generation X-ray telescope would 
be needed. We also note that this obj ect was detected in the 
analys is of the XMM-Newton data by Shtv kovskiv & Gilfanovl 
(2005), but not identified as a QSO. 

We also note that the narrow emission lines seen in the 
spectrum of QSO J050940.60-685320.6 near A rest = 3000Aare 
not intrinsic to the quasar. The line of sight to this object inter- 
sects an emission nebula and those features are spectral extrac- 
tion residuals from the [Om] lines. The strong emission line is 
real; lack of other emission features in the spectrum identifies 
it as Mg H/12798 at z em = 0.67. 

The analysis of the properties of these and other QSOs 
seen in X-rays will be presented in the forthcoming paper. The 
non-quasars among the X-ray-selected objects were mostly low 
mass X-ray binaries and hot stars. We will present the analysis 
of those and other X-ray objects observed during the course of 
our project in a separate paper. 

3. Properties of variability-selected quasars 

Combining objects presented here with objects identified in 
previous studies, we end up with fourteen quasars selected 
from the OGLE-II LMC and SMC database s based on their 
varia bility: in the SMC five QSOs identified bv lDobrzvcki et all 
(20033), and in the LMC: six identified in this paper, two 
identified by iGeha et aD ( l2003h and one by iDobrzvcki et alJ 
120021) . This sample is quite well defined: all those objects were 
selected using the same technique from homogeneous, good 
quality photometric data. We will present below the analysis of 
some observable properties of those objects. 

3.1. Colour-colour diagram 

Ii 1 I 
Dobrzvcki et al. (2003a) we found that variability-selected 

quasars could be distinguished from non-QSOs based on the 
V-I colour, since quasars were typically redder than other vari- 
able candidates. This effect is clearly seen with new quasars 
added. In Figure[3]we show the V-I vs. B-V colour-colour di- 
agram for all LMC and SMC candidates. One can clearly see 
that setting a cutoff at V-I ~ 0.3 would dramatically increase 
the efficiency of the variability-based selection method, to 40- 
50%. 

We note here that the thr ee objects w hich were not observed 
- L142, L143 and L146 in lEverl J2002I) - all have V-I < 0. 
Their location on Fig.|3]indicates that they are rather unlikely 
to be quasars. 

It is, of course, a very well known fact that quasars 
preferable occupy specific regions in colour-colour diagrams 



and this information was applied to 


quasar surveys (e.g. 


Cristiani et al. 


1989; Newberg & Yannvlll997l [Newberg et alJ 


1999; Richards et al 


. 2001). However, we stress that the di- 



agram presented in Fig. [3] only shows the distinction be- 
tween colour properties of stars and quasars which were 
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Fig. 3. The V-I vs. B-V colour-colour diagram for objects 
in the direction of the LMC and SMC. Filled symbols show 
QSOs. The preferred location for quasars can be clearly seen. 

QSO candidates, already preselected with colour-based crite- 
ria. Applicability of this plot is therefore limited to our sample. 

3.2. Brightness-colour relation 

Several papers (e.g.lPaltani & Courvoisierll994llHelfand et alJ 
120011: IVandenBerketalJl2004l and references therein) indi- 
cate existence of a brightness-colour relation in variable QSOs: 
quasars tend to be bluer in higher states. 

In addition to a few hundred measurements in the / band, 
the OGLE-II database provides a few tens of measurements in 
the B and V bands. It is therefore possible to investigate this 
effect in our sample. 

The diagram in Figure [4] shows the parameters of the / vs. 
V-I relation obtained from linear least square fits. The x-axis 
shows the Pearson's correlation coefficient while the y-axis 
shows the fitted slope and corresponding uncertainty. The ob- 
jects in which the "brighter-bluer" effect was significant would 
occupy the upper right corner of Fig.|4j and their slopes would 
be significantly greater than zero. 

One can clearly see that this effect is not generally present 
in our data. With the excepti on of two objects (both behind the 
SMC: SI and S12 in lEverll2002l) . most of our quasars do not 
show this effect and some even give an indication for opposite 
relation. This is somewhat curious since the effect is rather well 
established in the literature and one would expect that our set 
of objects was large enough to show it. 

We note that the behaviour of the B-V colour is quantita- 
tively similar. 

One cautionary note here is that sometimes the BVI spec- 
tral bands include QSO emission lines in our objects. This can 
weaken the effect, since the lines may account for a sizeable 



Dobrzycki et al.: Discovery of nine quasars behind the Large Magellanic Cloud 



5 




-1.0 -0.5 0.0 0.5 1.0 

Correlation coefficient 

Fig. 4. Parameters of the / vs. V-I relation: Pearson's corre- 
lation coefficient versus slope for all identified objects in the 
direction of LMC and SMC. Filled symbols show quasars. See 
text for discussion. 

fraction of the observed flux, while the effect is associated with 
continuum emission. Possible variability of emission lines - 
e.g. the Baldwin effect - may affect the result in our case. 

We also note that if this effect was indeed seen it could 
serve as an auxiliary selection criterion since variable stars 
clearly do not follow this relation. 

3.3. Amplitude-timescale relation 

Analysis of variability properties of the objects in the sample 
appears to show an interesting phenomenon. In Figure |5] we 
plot the residual variability dispersion in OGLE-II data versus 
the variogram slope (see Ever 2002) for all objects with known 
identifications. 

There is no correlation between those two quantities for 
QSO candidates that turned out to be stars. However, a de- 
cent correlation (R = 0.64) is seen for quasars, indicating 
that the magnitude of the QSO brightening is correlated with 
the timescale in which it occurs. In other words, unlike stars, 
quasars do not show large outbursts or long periods of relative 
calmness. 

We emphasize that the above statements only apply to stars 
in the QSO candidate sample and not to stars in general. A 
similar effect would be seen in several types of variable stars, 
which were removed from the original QSO candidate sample 
by magnitude and colour cuts. 

This effect should be verified on a larger sample. If it is 
confirmed, it can potentially be used to increase the efficiency 
of the candidate selection, by filtering out those of the stars that 
lie above and below the region in Figure[5]in which QSOs are 
found. 
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Fig. 5. Variability versus variogram slope for identified objects. 
Filled symbols show QSOs. 

4. Discussion 

The quasars presented in this paper have increased the number 
of known QSOs in the direction of LMC by 18%. On Figure^ 
we present the locations of the new quasars, together with 
the locations of previously known QSOs. In total, this picture 
shows 61 quasars. The image shows a large qualitative jump in 
the number of known objects since 2002. 

However, the area in the direction of the LMC bar appears 
to be underrepresented. This is somewhat surprising. The QSO 
candidates were distributed quite uniformly over the bar, while 
the newly discovered objects are all on its outskirts. It is pos- 
sible that some effects - for example, increased extinction or 
source crowding affecting variability analysis - do influence 
the candidate selection in the densest fields. However, since the 
number of object is still rather small, we are not in a position 
to determine whether a true bias is involved. 

We confirm the result seen in the SMC data 
JPobrzvcki et alJ l2003al) that the addition of the colour 
criterion to the variability selection procedure will lead to 
much higher efficiency in quasar searches, since QSOs are 
typically redder than other variability-selected objects. As 
expected, this distinction for the LMC objects was somewhat 
less strict compared to the SMC, but even here it would lead to 
very high, approaching 50%, success rates. 

In our data we see a marginally significant correlation be- 
tween the variability magnitude and the variability timescale. 
This correlation is entirely absent for variable stars from the 
QSO candidate pool. If this effect is confirmed, it can be used 
to increase the method efficiency even more. 

Quasars from our sample do not generally follow the 
"brighter-bluer" relation seen in larger samples, such as SDSS. 
However, this may be affected by variability of the QSO emis- 
sion lines, which enter observed photometry bands. 
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Fig. 6. Quasars in the vicinity of the Large Magellanic Cloud. The image is roughly 7° x 7°. Open squares and diamonds show 
the positions of, respectively, variability- and X-ray-selected quas ars pre sented in this paper. The filled squares show variability- 
selected QSOs from our candidate list confirmed by Geha et al. (2003). Th e "X" marks show the variability-selected quasars 
from lGeha et al.l d2003h . Open circles show four X-ray selected QSOs from Dobr zvcki et alJ J2002I) . Crosses show the quasars 
known before 2002; a dramatic improvement in both the number of known objects and the coverage is clearly seen. Only few 
QSOs are still known behind the Cloud's bar. LMC image courtesy of G. Bothun. 



We plan to utilise the variability selection method for 
searching for quasars behind the Galactic B ulge, in low extinc - 
tion windows monitored by OGLE-II (see ISumi e t al. 2005). 
A successful search would carry great hope for studies of the 
dynamics of the Bulge, since the expected proper motions are 
much larger than in the Magellanic Clouds. 
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